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Abstract 

We present calculations of the localisation length, A2, for two interacting par- 
ticles (TIP) in a one-dimensional random potential, presenting its dependence 
on disorder, interaction strength U and system size. \2(U) is computed by 
a decimation method from the decay of the Green function along the diag- 
onal of finite samples. Infinite sample size estimates ^(U) are obtained by 
finite-size scaling. For U = we reproduce approximately the well-known 
dependence of the one-particle localisation length on disorder while for finite 
U, we find that &(U) ~ 6(0) /3{C/) with (3(U) varying between (3(0) = 1 and 
(3(1) ~ 1.5. We test the validity of various other proposed fit functions and 
also study the problem of TIP in two different random potentials correspond- 
ing to interacting electron-hole pairs. As a check of our method and data, we 
also reproduce well-known results for the two-dimensional Anderson model 

without interaction. 
71.55.Jv, 72.15.Rn, 71.30.+h 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The interplay between disorder and many-body interactions in electronic systems has 
been studied intensively over the last two decades |l[ and still continues to receive much 
attention. Unlike the case of non-interacting electrons, where the "scaling hypothesis of 
localisation" can reliably predict the results of many experimental and numerical stud- 
ies, there is no equally successful approach to localisation when many-particle interactions 
become important M. Recently, experimental studies of persistent currents in mesoscopic 
rings and the discovery of a metal-insulator transition in certain two-dimensional (2D) elec- 
tron gases at zero magnetic field |§ have shown that the presence of interactions can indeed 
give rise to both quantitatively and qualitatively unexpected phenomena. 

A simple and tractable approach to the problem of interacting electrons in disordered 
materials is the case of only two interacting particles (TIP) in a random potential in one 
dimension. For a Hubbard on-site interaction this problem has recently attracted a lot of 
attention after Shepelyansky |4],|5j argued that attractive as well as repulsive interactions 
between the two particles (bosons or fermions) lead to the formation of particle pairs whose 
localisation length A2 is much larger than the single-particle (SP) localisation length Ai 
Based on a mapping of the TIP Hamiltonian onto an effective random matrix model (RMM) 
he predicted 

A 2 ~ {U/t) 2 \\ (1) 

at two-particle energy E = 0, with t the nearest-neighbor hopping matrix element and U 
the Hubbard interaction strength. Shortly afterwards, Imry || used a Thouless-type block- 
scaling picture (BSP) in support of this. The major prediction of this work is that in the 
limit of weak disorder a particle pair will travel much further than a SP. This should be 
contrasted with renormalization group studies of the ID Hubbard model at finite particle 
density which indicate that a repulsive onsite interaction leads to a strongly localised ground 
state 0. 
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The preferred numerical method for accurately computing localisation lengths in disor- 
dered quantum systems is the transfer matrix method (TMM) ||10|| . Thus it was natural that 
the first numerical studies devoted to the TIP problem also used the TMM to investigate 
the proposed enhancement of the pair localisation length A2 [PJIIfl. Other direct numerical 



approaches to the TIP problem have been based on the time evolution of wave packets f§,[[2 



exact diagonalization ||13| , or Green function approaches |jT4| , |l5|l . In these investigations usu 



ally an enhancement of A2 compared to Ai has been found but the quantitative results tend 
to differ both from the analytical prediction in Eq. ([!]), and from each other. Furthermore, a 
check of the functional dependence of A2 on Ai is numerically very expensive since it requires 
very large system sizes. 



Following the approach of Ref. [11 1, two of us studied the TIP problem by a different 



TMM ]rjj at large system size M and found that (i) the enhancement A2/A1 decreases with 
increasing M, (ii) the behavior of A2 for U = is equal to Ai in the limit M — > 00 only, and 
(iii) for U ^ the enhancement A2/A1 also vanishes completely in this limit. Therefore we 
concluded jlBJ that the TMM applied to the TIP problem in ID measures an enhancement 
of the localisation length which is due to the finiteness of the systems considered. The main 
problem with this approach is that the enhanced localisation length A2 is expected to appear 
along the diagonal sites of the TIP Hamiltonian, whereas the TMMs of Refs. [|ll|)IH] proceed 
along a SP coordinate. Various new TMM approaches have been developed to take this into 



account [|IT| , |16| -|T8|| , but still all TMMs share a common problem: in general the U = result 
for A2 does not equal the value of X\/2 which is expected for non-interacting particles as 
explained below. Rather, they show localisation lengths \2(U = 0) which are much larger 
than Ai/2 and very close to \2(U ~ 1). 

The obvious failure of the TMM approach to the TIP problem in a random potential 
has lead us to search for and apply another well-tested method of computing localisation 
lengths for disordered system: the decimation method (DM) |T3. Furthermore, instead of 



simply considering localization lengths A 2 (t/) obtained for finite systems [ pH| , pT^ T5| , or by 



simple extrapolations to large M | lq| , we will construct finite-size scaling (FSS) curves and 



compute from these curves scaling parameters which are the infinite-sample-size estimates 
of the localization lengths ^(U). We find that onsite interaction indeed leads to a TIP 
localisation length which is larger than the SP localisation length at E = and for not too 
large U. However, the actual functional dependence is not simply given by Eq. ([I]). In fact 
our data allow us to see ^(U) ~ £2(0)^ with an exponent (3 which increases with increasing 
\U\ at E = 0. 

The paper is organized as follows: In section |T| we introduce the numerical DM used 
to compute the localisation lengths. In section |T|, we investigate the numerical reliability 
of the DM by studying the Anderson model in 2D. We then apply the method to the case 
of TIP in section [IV] and use FSS in order to construct infinite-sample-size estimates in 
section |V|. We fit our data with various functional forms for £ 2 put forward in the literature. 
In section [VI] we also study the problem of two interacting particles in different random 
potentials. In section [VI H we study the related problem of a SP in a 2D random potential 
with additional barriers. We conclude in section |V111| . 

II. THE DECIMATION METHOD 

We shall be considering properties of Hamiltonians of the form 

H = — t (\n,m)(n + l,m\ + \n,m)(n,m + 1| + h.c.) 

n,m 

+ K m H e n + e « + U(n)6 nm ) (n,m\ (2) 

where the choice of e*, and the definition of U(n) depends on the specific problem 
considered. For the case of TIP in ID the indices n and m correspond to the positions of 
each particle on a ID chain of length M and e\ = e\ G [-W/2, W/2]. We shall also present 
results for the case of 7^ which corresponds to two interacting particles in different ID 
random potentials, e.g., two electrons on neighboring chains, or an electron and a hole on the 
same chain. In these cases U(n) = U is the interaction between the two particles. Instead of 
considering TIP we can also choose M 2 uncorrelated random numbers e nm £ [—W/2, W/2] 
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and replace + in (||) by e nm . Then the Hamiltonian (0) corresponds to the standard 
Anderson model for a single particle in 2D with an additional potential U(n) along the 
diagonal of the 2D square. In all cases we use hard-wall boundary conditions and t = 1 sets 
the energy scale. 

We now proceed to construct an effective Hamiltonian along the diagonal of the M x M 
lattice by using the DM If we write A(E) = El— H, the defining equation A(E)G(E) = 



1 for the Green function G(E) can be written as 

jV-l 

M E )G 3 k(E) + A lN {E)G Nk {E) = 5 lk (3) 

5=1 

where N = M 2 is the total number of sites in the system and the indices i, j, k = 1, . . . , N 
represent multi-indices for the M 2 states \n,m). From this we can see by choosing i = N 



that 



N-l 



g »( b > = j%-£#t§ g ^>- (4) 

Ann[E) ~[ A NN {ti) 



Substituting into (0) gives for all k ^ N 

E W E ) - Am{ ;>'Tjr' Gjk(E) = 5 lk . (5) 

5=1 L 



JV ' A NN (E) 

In this way we have obtained an effective Hamiltonian H'(E) with matrix elements H[-{E) = 
Hi a + ^inH N j w j lose Q reen function is identical to that of the full Hamiltonian on all non- 

l E-H NN 

decimated sites. This process is repeated until we are left with an effective Hamiltonian for 
the doubly-occupied sites only. We remark that due to cpu-time considerations it turned 
out to be useful to split the Hamiltonian into two halves along the diagonal and to start the 
decimation process from the outer corner of the triangular half and then decimate in slices 
towards the diagonal. The procedure is shown pictorially in Fig. [I]. Furthermore, for the 
case of TIP we only need to decimate one half and can use the symmetry of the spatial part 
of the wave function for the other half. 

We shall now focus our attention upon the TIP localisation length A2 obtained from the 
decay of the transmission probability of TIP from one end of the system to the other. In 



accordance with the SP case QTO | , A2 is defined by the TIP Green function, G%(E). More 
precisely ]14[] 

^ = -|M3Tj ln l<M|G 2 |M,M>|. (6) 

The Green function matrix elements (1, l|Ga|M, M) are computed by inverting the matrix 
A(E) = El — H(E) obtained from the effective Hamiltonian 11(E) for the doubly occupied 
sites. We remark that in order to reduce possible boundary effects, we compute A2 by 
considering the decay between sites slightly inside the sample instead of the boundary sites 
(M, M). 

III. TESTING THE DECIMATION METHOD 

As mentioned in the introduction, one of the surprises of the TIP problem is the apparent 
inapplicability of the TMM approach, which leads to large enhancement of the localisation 
lengths even in the absence of interaction (U = 0). Thus it appears necessary that before 
using the DM for the case of TIP, we should also check that by restricting ourselves to 
the decay of the Green function along the diagonal, we do not encounter similar artificial 
enhancements of \%(U = 0). As a first test, we have therefore studied the decay of the 
Green function along the diagonal for the usual 2D Anderson model at various disorders 
W = 0.65, . . . , 20 and system sizes M = 51, . . . , 261. For comparison, estimates of Ai were 



computed by the standard TMM |10] in 2D. We then use FSS as in |I(J and compute the 
localisation lengths £i(W) valid at infinite system size for both sets of data. 

In Fig. we show the resulting localisation lengths £1 obtained by TMM with 1% ac- 
curacy and DM averaged over 100 samples for each W and M. When we are considering 
a 2D system, to obtain the correct value of the localisation length we have to multiply the 
localisation length obtained from Eq. (|6]) by y/2 to take account of the fact that we are 
studying the decay along the diagonal. As shown in Fig. 0, the agreement is good down to 
W = 4.5 where the FSS becomes unreliable. We clearly see that using the DM to calculate 
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the Green function along the diagonal reproduces the well-known results of Ref. |H| up 
to a geometrical factor which is easily understood. Furthermore, the deviations from the 
TMM results for £1 show that our method underestimates the infinite system size results. 
Therefore the above mentioned problem of the TMM giving rise to too large a value for 
the TIP-localisation lengths A2 due to small system size should not appear. We emphasize 
that the FSS procedure is more than an extrapolation to the infinite system size |10[ and it 
allows us to identify the disorders at which FSS breaks down as shown in Fig. ||]. 

Before proceeding to the case of TIP, we need to discuss an important difference between 
the data obtained from TMM and DM. The TMM proceeds by multiplying transfer matrices 
for 2D strips (3D bars) of finite size M (M x M) many times until convergence is achieved. 
The localisation lengths are then computed as eigenvalues of the resulting product matrix 
||10|| . However, in the present case of DM (or any other Green function method applied to 
TIP), the localisation lengths are estimated by assuming an exponential decay as in Eq. 

Such a simple functional form, however, will no longer be reliably observable when 
£1 ~ M and we will start to measure the oscillations in the Green function underlying the 
exponential envelope (^). Looking at Fig. Q we indeed see that the deviations from the 
TMM result start at £1 ~ 250, that is, just at the largest system sizes used. Increasing the 
number of samples will reduce this effect, but this quickly becomes prohibitive due to the 
immense computational effort. With this in mind, we now continue to the case of TIP. 



IV. THE TIP PROBLEM AT FIXED SYSTEM SIZE 

We now compute the Green function at E = for 26 disorder values W between 0.5 
and 9 indicated in Fig. ^|, for 24 system sizes M between 51 and 251, and 11 interactions 
strengths U = 0, 0.1, ... , 1.0. For each such triplet of parameters (W, M, U) we average the 
inverse localisation lengths I/A2 computed from the Green function according to Eq. 
over 1000 samples. 

In Fig. [| we show the results for M = 201. Let us first turn our attention to the case 
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U = 0. As pointed out previously |15], the TIP Green function G2 at E = is given by 



a 



convolution of two SP Green functions Gi at energies Ex and —E 1 as 

(l,l|G a (0)|M,M)~ J dE'(l\Gx{E')\M){l\Gx{-E')\M). (7) 

Assuming that (l\Gi(E)\M) oc exp[-|M- l\/\x(E)], where X X (E) is the SP localisation 
length of states in the ID Anderson model 0, one expects that the largest localisation 
lengths dominate the integral. Since Ai(0) > Xi(E), this implies that (1, 1|G 2 (0)|M, M) « 
exp [— 2|M — 1 1 / Ai (0)] . Applying Eq. @, we get A2 = Ai/2 p0| . Therefore we have also 
included data for Ai/2 in Fig. |3|. Since Ai deviates from the simple power-law prediction |7j 
Ai « 104/ at E = for A x < 4 (W > 5), we have computed A x by TMM in ID with 
0.1% accuracy. 

Comparing these results with the TIP localisation lengths obtained from the DM, we 
find that for 1 < W < 6, the agreement between \2{U = 0) and Ai/2 is rather good and, 
contrary to the TMM results, there is no large artificial enhancement at U = 0. For smaller 
disorders W < 1, we have A2 ~ M/2 so that it is not surprising that the Green function 
becomes altered due to the finiteness of the chains |2T| . This results in reduced values of A 2 . 
For large disorders W > 6, we see a slight upward shift of the computed A 2 values compared 
to Ai/2. This effect is due to a numerical problem, since straightforward application of Eq. 
(§) is numerically unreliable for values of A 1 as small as 1. 

It is noticeable from these results, however, that the values of \2(U = 0) are still slightly 
larger than Ai/2. In order to explain this behavior, we have computed (l\Gi(E)\M) by exact 
diagonalization of the SP Hamiltonian for at least 100 samples at many different energies 
inside the band and then integrated as in Eq. (0). Plotting the resulting localisation lengths 
in Fig. R] we see that indeed the agreement with A 2 (t/ = 0) is better than with Ai/2. Thus 



the corresponding conjecture of Ref. [|T^] is shown to be true. 

For U between 0.1 and 1 we have found that the localisation lengths are increased by the 
onsite interaction (cp. Fig. We have also seen that for W > 1.4 the localisation lengths 
A 2 (t/) increase with increasing U. For smaller W we have A 2 (C7) ~ M/2 and, as discussed 
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above, the data become unreliable for fixed system size. 

Up to now we have been mostly concerned with the behavior of A2 as function of disorder 
for U G [0, 1]. However, for large U, it is well-known that the interaction splits the single TIP 
band into upper and lower Hubbard bands. Thus we expect that for large U the enhancement 
of the TIP localisation length vanishes. In Fig. f|we present data for A 2 (£/)/A 2 (0) obtained 
for three different disorders for system sizes M = 201 at E = 0. In agreement with the 



previous arguments and calculations jl6| , |22| , |23| , we find that the enhancement is symmetric 
in U and decreases for large \U\. For small \U\, we see that the localisation length increases 
nearly linearly in \U\ with a slope that is larger for smaller W. We do not find any U 2 
behavior as in Refs. |§|j]|§. In Ref. p3) is has been argued that at least for Ai ~ M, 
there exists a critical U c = 24 1//4 , which is independent of W , at which the enhancement 
is maximal. We find that in the present case with Ai < M the maximum enhancement 
maxj/ [A 2 (£7)/A 2 (0)] depends on the specific value of disorder used. Another observation of 
Ref. is the duality in U and y/24/U for very large \U\ (small 1/\U\). The data in Fig. | 
are only compatible with this duality for the large disorder W = 5. For the smaller disorders 
and for the range of interactions shown, we do not observe the duality. We emphasize that 
this may be due to restricting ourselves to values U < 4. 

For E 7^ 0, the independence of the enhancement on the sign of the interaction U is 
no longer valid. In Fig. |^ we show A^LQ/A^O) for the same disorders as before, but now 
at energies E — ±1. We find that the enhancement for U = 1 is larger at E — 1 than at 
E = — 1, whereas exactly the opposite is true for U = — 1. Thus we see that for positive 
(negative) U the energies of TIP states are shifted towards higher (lower) values, eventually 
leading to the formation of the aforementioned Hubbard bands. In Fig. || we show the 
localisation lengths at several values of E for W = 4. As expected from the discussion above 
the localisation lengths are always smaller than at the band center. The enhancements, 
however, which are shown in Fig. |5|, can be equally large for E — and E ^ 0. 
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V. FSS APPLIED TO THE TIP PROBLEM 



In order to overcome the problems with the finite chain lengths, we now proceed to use 
the FSS technique and construct FSS curves for each U — 0, 0.1, . . . , 1. In Fig. [7] we show 
the data for the reduced localisation lengths A2/M which is to be rescaled just as in the 



standard TMM to obtain the localisation length £ 2 for the infinite system. Note that 
data for small W is rather noisy and will thus most likely not give very accurate scaling. 
Furthermore, in Fig. |8| we show A2 for W = 3 and W = 9 and all interaction strengths 
{7 = 0, 0.1, . . . 1.0. We see that whereas for W = 3 the values of A2 for {7 = show only 
small variations for large M, the W = 9 data shows a rapid increase of A2 as M increases. 
This is due to the numerical problem of estimating a small localisation length of the order 
of 1 in a large system by Eq. (||). It is most pronounced for small U where the localisation 
lengths are the smallest. Going back to Fig. [7], we see that this does not influence the 
reduced localisation lengths A 2 /M very much and thus is not expected to deteriorate the 
FSS procedure. However, in order to set an absolute scale in the FSS procedure, one usually 
fits the smallest localisation lengths of the largest systems to A2/M = ^/M + b(C,2/M) 2 with 
b small [HJ. In the present case this would mean taking the unreliable data for W = 9. 
Therefore, for each U we fit to the localisation length at W = 3 and adjust the absolute scale 
of £ 2 accordingly. In Fig. [] we show the resulting scaling curves A 2 /M = f^/M) for {7 = 0, 
0.2 and 1.0. Note that, as expected from Fig. 0, FSS is not very accurate for small W. 
The previously discussed unreliable data for large W are visible only in very small upward 
deviations from the expected 1/M behavior. In Fig. ^ we show the scaling parameters £2 
obtained from the FSS curves of Fig. ^. 

A simple power-law fit £2 oc W~ 2a in the disorder range W G [1,5] yields an exponent a 
which increases with increasing U as shown in the inset of Fig. [10], e.g., a = 1.1 for U = 
and a = 1.55 for {7=1. Thus, although in Fig. |3| the A2 data at M — 200 nicely follows X±/2 
for U = 0, we nevertheless find that after FSS with data from all system sizes, £2(0) still 
gives a slight enhancement. Because of this in the following we will compare £2 ({7 > 0) with 
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£2(0) when trying to identify an enhancement of the localisation lengths due to interaction. 
We emphasize that the slight dip in the a{U) curve around U = 0.7 has also been observed 
in Ref . . 



The derivation of Eq. ([[]) is based on a mapping of the TIP Hamiltonian onto an effective 
random matrix model while assuming uncorrelated interaction matrix elements M. In Refs. 



21 and 221 a more accurate estimate of the matrix elements of the interaction in the basis of 



SP states was calculated showing that the original estimates of Ref. Q were oversimplified. 



The authors of Ref. [£2| then considered a more appropriate effective random matrix model 
and obtained A2 oc \f for large values of Ai. To correct for smaller values of Ai they 
suggested a more accurate expression should be A 2 oc Af (1 + c/Ai). An important prediction 
of this work is that j3 is [/-dependent with f3 ranging from 1 at small U and very large U 
to nearly 2 for intermediate values U oc t. Using our data obtained from FSS, we translate 
this fit function into 

&V) oc 6(0)^1 + -^rV (8) 



6(0)7 ' 

We remark that the actual least-squares fit is performed with the numerically more stable 
fit function y = a + (3x + c * exp(— x) with y = ]n[£ 2 (U)] and x = ln[£ 2 (0)]. In Fig. |TT] we 
show results for disorders W G [1, 6] and various U. As can be seen easily, the fit is rather 
good and does indeed capture the deviations from a simple power-law 6(^) oc £2(0/ for 
small localisation lengths. In the inset of Fig. [11] we show the variations of j3 with U for 
both the simple power-law and the fit according to Eq. (]8|). We note that contrary to Ref. 



22], we find (3 < 1.5 for all U values considered. 



In Ref. |TJ] is has been suggested that a more suitable functional dependence of the TIP 
localisation lengths is given by A2 = Ai/2 + c|C/|Af. Using the £2 data and taking instead of 
Ai/2 the more suitable £2(0) we translate this proposed fit as 

UU) - 6(0) oc 6(0)^. (9) 

In Fig. we plot £ 2 (U) - f 2 (0) vs. f 2 (0) for [/-values 0.2 and 1. We find that instead of 
being able to fit the data with a single (5, it appears that for small 6(0) < 10 we have f3 ~ 2, 
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whereas for larger £2(0) we find (3 ~ 3/2. Note that a crossover from the functional form 
(0) with (3 = 2 to f3 = 3/2 has been suggested previously However, in that work, the 
exponent 3/2 is supposed to be relevant for larger disorders, opposite to what we see here. 
As pointed out previously, our FSS may give rise to artificially small values of £z(U) close 
to the largest system size, and one might want to argue that the reduction in slope is due to 



this effect. However, we emphasize that the crossover observed in Fig. [12] occurs at W = 2.5 
where FSS appears to be still reliable. We remark that an exponent close to 1.5 for small 
W has also been found in Ref. [25| from a multifractal analysis. 



The most recent suggestion of how to describe the TIP localisation data is due to Song 
and Kim fL5| . They assume a scaling form 

6 = w-e° g (\u\/w A ) (10) 

with g a scaling function and obtain A = 4 by fitting the data. Choosing the same value for 
A we find that our data can be best described when (3q is related to the disorder dependence 
of £2 as (j3 — /3 )/A 1/4. However, the scaling is only good for W G [1,5] and U > 0.3. 
Unfortunately, even using our varying exponent /3{U), we have not been able to obtain a 
good fit to the scaling function with the data for all U. We emphasize that the £2 values for 
U < 0.2 are smaller than for U > 0.3 and thus numerically quite reliable. 
A much better scaling can be obtained when plotting 

UU) - 6(0) = g [/(l%(0)] (11) 



with f(U) determined by FSS. In Fig. [13] we show the resulting scaling curves and scaling 
parameters f(U). Note that the scaling is valid for U = 0.1,0.2,.. .1.0 and most disorders 
W G [0.6,9]. Again we see the crossover from a slope 2 to a slope 3/2. Deviations from 
scaling occur for large and very small values of &(U) and are most likely due to numerical 
inaccuracy as discussed before. The behavior of f(U) as shown in the inset indicates that 
for U > 0.3 a linear behavior f(U) oc U may be valid which then translates into U 2 (U 3 ^ 2 ) 



dependence of ^(U) — £2(0) in the regions of Fig. [13| with slope 2 (3/2). However, for 
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U < 0.5, one could also argue that f(U)ocVU which would give £ 2 (C/) - f 2 (0) oc U (f/ 3/4 ) 
in these regions. We note that a crossover from U to U 2 behavior had been proposed in Ref. 



24J, but it should appear at larger values of U and also be W dependent. We observe that 
the best fit to the f(U) data is obtained by a logarithmic ^/-dependence as indicated in the 
inset. 

Thus in summary it appears that our data cannot be described by a simple power-law 
behavior with a single exponent neither as function of W, nor as function of £2(0), nor after 
scaling the data onto a single scaling curve. The best power-law fit is obtained in Fig. |ll| 
with an exponent /3(U), whereas after scaling of ^(U) — £2(0) onto a single curve we need 



at least two powers to describe the scaling curve as shown in Fig. [L3| Lacking a convincing 
explanation as to what fit function should be correct, we must at present be content with 
letting the reader decide for himself. 

VI. THE INTERACTING ELECTRON-HOLE PROBLEM 

Let us now consider what happens when the two particles are in different random poten- 
tials such that in general e\ 7^ e^. Such a problem is relevant for the proposed experimental 
verification of the TIP effect by optical experiments in semiconductors |26[. In these ex- 
periments, the electron will be in a random potential different from that of the hole. Thus 
this choice of random potential models the case of interacting electron-hole pairs (IEH). 
Again, we will mostly be concerned with the case of repulsive interactions. In the experi- 



mental situation, of course, the interaction is attractive. As shown in Fig. [L4] we again have 
\z(—U) = A 2 (6 r ) for E = and thus our results apply also to the case U < 0. For simplicity, 
we also take the width of the disorder distribution to be the same for both particles. 

As for TIP we compute the IEH localisation lengths by the DM along the diagonal. 
Comparing with the results presented in the previous sections, we find that the results for 
IEH are very similar to the case of TIP. FSS is possible and again the best fit is obtained by 
using Eq. (H) as shown in Fig. O for U = 0, 0.1, 0.2, 0.5 and 1.0. The values of the power 
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/3{U) shown in the inset of Fig. [15] are also much as before. Thus we can conclude that the 
case of IEH is very close to the TIP problem. 



VII. THE 2D ANDERSON MODEL WITH AN ADDITIONAL DIAGONAL 

POTENTIAL 

In Ref. two of us argued that straightforward application of the random matrix 
models (RMM) J| and the block-scaling picture (BSP) || gives rise to an erroneous en- 
hancement of the SP localisation length ^ in a 2D Anderson model with additional random 
perturbing potential U(n) G [-U, U] along the diagonal. In fact, the same is true if the 
potential along the diagonal is taken to be constant, i.e. U(n) = U. Although it appears 
obvious that no such SP enhancement should exist, we have checked it here with the DM. 
In Fig. [16| we show examples of the resulting SP localisation lengths £1 obtained as before 
from FSS of SP localisation lengths Ai calculated for various system sizes M = 51, . . . , 261, 
disorders W and potentials U = 0, 0.1, . . . , 1. As expected, we find that for large disorders 
W > 5, the data is well described by the 2D TMM results already presented in section fT\\ . 
There are only small changes due to the additional random potential, all of which tend to 
decrease the localisation lengths as they should. This is in contrast to the straightforward 
application of the RMM and the BSP |yj which therefore fail for the 2D SP Anderson model 
with additional random potential along the diagonal. Of course this does not mean that 
these methods also have to fail for TIP, where, as we have shown in the previous sections, 
a tendency towards delocalisation due to interaction definitely exists. 



VIII. CONCLUSIONS 

In conclusion, we have presented detailed results for the localization lengths of pair 
states of two interacting particles in ID random potentials. By using the DM to calculate 
the Green function along the diagonal it is possible to consider the 2D Anderson model and 
the problem of two interacting particles in ID within the same numerical formalism. We 
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have checked that for the 2D Anderson model without interaction the infinite system size 
results obtained via FSS from the DM data are in good agreement with results obtained 
from the standard TMM especially for localisation lengths up to the largest system sizes we 
have considered. It is also apparent that the DM data deviate from the TMM only towards 
smaller localisation lengths and hence no artificial enhancement of localisation lengths due 
to the DM approach is expected. 

For TIP in ID we observe an enhancement of the two-particle localisation length up to 
75% due to onsite interaction. This enhancement persists, unlike for TMM, in the limit of 
large system size and after constructing infinite-sample-size estimates from the FSS curves. 
We have tried to fit our results to various suggested models. The best fit was obtained with 
Eq. (H) in which the enhancement ^(^O/^O) depends on an exponent (3 which is a function 
of the interaction strength U . Such a [/-dependent exponent had been previously predicted 



in Ref. p2| for interaction strengths up to U = 1 with (3 up to 2. However, we find that (3 
reaches at most 1.5 for U — 1. Thus we do not see a behavior as in Eq. ([I]) with exponent 
2 when using the fit function of Ref. [2~2~|. On the other hand, after scaling the data onto a 
single scaling curve and using the fit function (|9|) as proposed with (3 = 2 in Ref. |TJ 



we 

find indeed (3 = 2 for not too small disorder strength, e.g., W > 2.5 for U = 1, but observe 
a crossover to a behavior with (3 = 3/2 for smaller W. For values of U > 1.5 we observe 
that the enhancement decreases again; the position of the maximum depends upon W. Very 
similar results are produced by placing the two particles in different potentials which is of 



relevance for a proposed experimental test of the TIP effect [26 



As a final check on our results we consider the effect of an additional on-site potential 
(both random and uniform) on the results for the SP 2D Anderson model. As one may 
expect for the case of an additional random potential one observes only a small decrease in 
the localisation length while for an additional uniform potential there is a small change in 
£i towards decreasing localization lengths for positive U. 
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(c) (d) 




FIG. 1. The decimation process: (a) The lattice is split into two parts, (b) Each half is then 
'decimated' independently. (c,d) Sites (•) and nearest-neighbor hops (— ) in the original lattice are 
replaced successively by effective long-range hops (— • •— ) between the effective sites (0). (e) This 
decimation continues until (f) the diagonal is reached, (g) Finally, the two halfs are recombined. 
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FIG. 2. Small symbols (with left and bottom axis) denote the FSS curve used to compute the 
£i values. Large symbols (with right and top axis) indicate SP localisation lengths y/2 £i obtained 
by DM along the diagonal (o) and £i computed by TMM of quasi-lD strips (dashed line) after 
FSS. The filled symbols correspond to a disorder at which FSS appears to be unreliable. 
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FIG. 3. Two-particle localisation length A2 at energy E = for system size M = 201 and 
interaction strength U = (O); U = 0-2 (^) and U = 1 (□). The thick solid line represents ID 
TMM data for SP localisation length Ax/2, the dashed line is computed from the convolution of 
SP Green functions in Eq. (0). The thin line is the perturbative result Ai/2 ~ 52 /W 2 . 
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FIG. 4. Enhancement \2(U) / \2(0) for TIP as a function of interaction strength U at E = 
for disorder W = 3 (+), W = 4 (*), and = 5 (x) and M = 201. The data are averaged over 
100 samples. The thick (thin) lines indicate data for U > (U < 0), full (open) circles denote the 
maximum for each disorder. 
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FIG. 5. Enhancement \2(U) / ^(0) for TIP as a function of interaction strength U at E = ±1 
for disorder IP = 3 (+), VP = 4 (*), and VP = 5 (x) and M = 201. The thick (thin) lines indicate 
data for E = 1 (E = -1). 
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FIG. 6. TIP localisation length A2 as a function of interaction strength U at W = 4 for 
E = -2 (<), E = —1 (V), E 1 = (*), E 1 = 1 (A), and E = 2 (>) and M = 201. The thick (thin) 
lines indicate data for E 1 > (E < 0). 
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FIG. 7. Reduced TIP localisation lengths A 2 /M for U = (left) and U = 1 (right) for all 
disorders W and system sizes M obtained by averaging 1000 samples for each triple (U, W,M). 
Different letters indicate different disorders. 
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FIG. 8. TIP localisation lengths A 2 for W = 3 (left) and W = 9 (right) for U = 0, 0.1, . . . , 1 
from bottom to top. We remark that we have taken the same set of random numbers for all U to 
increase the numerical efficiency. This is probably the reason why for different U the fluctuations 
in the dependence of A2 on M are similar. 
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FIG. 9. Finite-size scaling plot of the reduced TIP localisation lengths \ijM for U = (O)i 
U = 0.2 (O) and U = 1 (□). The data for U = 0.2 (U = 1) have been divided by 2 (4) for clarity. 
Data corresponding to W = 1 are indicated by filled symbols. The two curves at the bottom show 
the data for U = 0.2 and 1 and W < 2.5, shifted downward by one order of magnitude for clarity, 
but with the data W < 1 fitted with scaling parameters obtained from the fit in Fig. [ll]. 
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FIG. 10. TIP localisation lengths 6 after FSS for U = (O), U = 0.2 (O) and U = 1 
(□). The solid line represents ID TMM data for SP localisation lengths Ai/2, the dashed lines 
indicate power-law fits. Inset: Exponent a obtained by the fit of £2 on W~ 2a to the data for 
U = 0,0.1,...,!. 
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FIG. 11. TIP localisation lengths &(U) after FSS for U = (O), U = 0.1 (A), U = 0.2 
(O), U = 0.5 (*) and U = 1 (□) plotted versus f 2 (0). The data are for W G [1,6]. The dashed 
lines show fits according to Eq. (||), the solid line sets the reference for [7 = 0. Inset: Exponent (3 
obtained by the fit of Eq. (||) to the data for U = 0,0.1, . . . , 1. The open symbols correspond to 
the fit with c = 0. 
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FIG. 12. TIP localisation lengths £ 2 plotted according to Eq. (g) for U = 0.2 (O) and (7 = 1 
(□). The solid line indicates a slope 2, the dashed line a slope 3/2. The filled symbols correspond 
to W = 1. 
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FIG. 13. Scaling plot of Eq. (£]]) with TIP localisation lengths £ 2 for all U and W G [0.6,9]. 
The solid line indicates a slope 2, the dashed line a slope 3/2. Inset: The values of f(U) needed 
to make the data collapse onto the U = 0.1 curve. Solid, dashed and dot-dashed lines are fits of 
f(U) for U > 0.3, U < 0.5, and all U, respectively. 
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FIG. 14. Enhancement X 2 (U)/X 2 (0) for IEH 
for disorder W = 3 (+), W = 4 (*), and W = 5 
100 samples. The thick (thin) lines indicate data 
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as a function of interaction strength U at E = 
(x) and M = 201. The data are averaged over 
for U > (U < 0). 
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FIG. 15. IEH localisation lengths &(U) after FSS for U = (O), U = 0.1 (A), U = 0.2 (O), 
U = 0.5 (*) and U = 1 (□) plotted versus £2(0). The lines are fits as in Fig. [ll|. Inset: Exponent 
(3 obtained by fitting Eq. (§) to the data for U = 0, 0.1, . . . , 1. 
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FIG. 16. SP localisation lengths as in Fig. | for U = (- - O) an d Ci obtained by DM 
for U = 1 with additional random (□) or constant (x) potential energies along the diagonal. 
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